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Abstract Convincing evidence of a past interaction between two rich clusters 
A399 and A401 was recently found by the X-ray imaging observations. In this 
paper we examine the structure and dynamics of this pair of galaxy clusters. A 
mixture-modeling algorithm has been applied to obtain a robust partition into 
two clusters, which allows us to discuss the virial mass and velocity distribution 
for each cluster. Assuming that these two clusters follow a linear orbit and they 
have once experienced a close encounter, we model the binary cluster as a two- 
body system. As a result, four gravitationally bound solutions are obtained. The 
recent X-ray observations seem to favor a scenario in which the two clusters with 
a true separation of 5.4ft," 1 Mpc are currently expanding at 583 km s _1 along the 
direction with a projection angle of 67°. 5, and they will reach a maximum extent 
of 5.65 h- 1 Mpc in about 1.0 /i" 1 Gyr. 
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1 INTRODUCTION 

According to the hierarchical bottom-up scenario, clusters of galaxies are thought to be formed 
by accreting and merging of subunits. The structure and dynamics of the rich galaxy clusters 
with ongoing merger events are of great importance for understanding the cluster evolution. 
N-body numeric simulations show that substructure may occur in individual rich clusters before 
their final collapse and virialization (White 1976; Burns et al. 1994). Since the cluster merging 
events are the most common and energetic phenomena in the Universe, more and more obser- 
vational efforts in optical and X-ray bands have been devoted to the nearby rich clusters with 
significant substructures (e.g., White, Briel & Henry 1993; Gastaldello et al. 2003). 
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When these subsystems are slightly more separated they may be classified as separate galaxy 
clusters. The interacting system of clusters A399 and A40f is a good example. They have long 
been treated as a merging pair of clusters since they are close to each other, with a projected 
separation of ~ 0.6° between their central cD galaxies (McGlynn & Fabian 1984; Oegerle & Hill 
1994). The temperature map for this binary cluster, derived from the ASCA spatially resolved 
spectroscopic data, possibly suggests a physical link or a massive dark matter filament between 
these two clusters (Markevitch et al. 1998). The X-ray observation with the ROSAT High 
Resolution Imager (HRI) unveiled a significant linear structure of A399 which points directly 
to the core of A401, and this feature might be resulted from a past violent interaction (Fabian, 
Peres & White 1997). Recent XMM-Newton observations also confirmed the enhanced X-ray 
flux and temperature in the region between two clusters (Sakelliou & Ponman 2004). 

Therefore it is of great interest to search for the optical anomalies in dynamics of the 
member galaxies. In general, the merger history can be modelled on the basis of the structure 
and dynamics of cluster galaxies, intracluster gas, and intergalactic medium (e.g., Colless & 
Dunn 1996). This paper will use the existing redshift measurements to investigate the possible 
structures in 2-dimension map and in velocity space. A prevalent mixture-modeling algorithm, 
known as the KMM algorithm (Ashman et al. 1994), will be applied to obtain a robust separation 
of these two clusters. Disregarding the rotation of the system, we will try to model this cluster 
pair as a two-body system on the basis of the velocity distribution and virial mass estimates. 
In §2, we present the spatial distribution and localized variations in velocity distribution for 
all the member galaxies in the A399/A401 system as a whole. We apply the KMM partition 
algorithm and discuss the velocity distribution for each cluster in §3. Then, the virial mass 
estimate and two-body modeling for this binary system of galaxy clusters are performed in §4. 
Finally, a summary is given in §5. 

2 PROPERTIES OF THE SAMPLE 

A399 and A401 are nearby clusters of galaxies (z ~ 0.07181 and 0.07366; Oegerle & Hill 2001), 
classified respectively as type I-II and I clusters by Bautz & Morgan (1970). With respect to 
the center of this binary system (2 h 58 m 30 s , +13°20"; J2000.0), 1331 extragalactic objects with 
positional offsets less than 100.0 arcmin were extracted from the NASA/IPAC Extragalactic 
Database (NED). However, only 240 galaxies have spectroscopic redshifts. The remainder ap- 
pear only in one of the imaging surveys from radio, infrared, and X-ray wavebands. 

Most of the spectroscopic data were contributed by Hill & Oegerle (1993) who carried out 
a survey of the cD clusters. The redshift measurements were detailed in Hill & Oegerle (1993), 
and the typical velocity uncertainty for the galaxies is less than 100 km s^ 1 . The distribution 
of spectroscopic redshifts for 240 known galaxies is shown in Fig. 1. There are 215 galaxies 
with their redshifts in a range of 18,000 km s~ 1 < cz < 25,000 km s _1 , with a peak at ~ 
21,500 km s _1 . Only one peak can be found in the velocity distribution because the velocity 
dispersions for individual clusters are larger than the apparent velocity separation between 
two clusters. It is unambiguous to treat these 215 galaxies as the members of this cluster pair 
since their redshift distribution is quite concentrated and isolated. The contamination from the 
foreground and background galaxies should be very slight. The spatial distribution for these 
215 galaxies is presented in Fig. 2. We superpose the contour map of the surface density that 
has been smoothed by a Gaussian window with a = 2'. Because of the severe overlap in the 
redshift distributions, clusters A399 and A401 can not be completely resolved by the velocity 
distribution only. 

In order to detect the substructures in both the velocity space and the projected map, 
we take use of the K-test (Colless & Dunn 1996) for the A399/A401 system as a whole. The 
statistic n n is defined to characterize the local deviation on the scale of groups of n nearest 
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Table 1 Result of K-Test for 215 Member Galaxies in the Binary System 



n 


4 


5 


6 


7 


8 9 


10 


11 


12 


P^Kri ^ f^n ) 


18.3% 


44.4% 


42.6% 


34.1% 


17.0% 11.9% 


6.8% 


15.1% 


14.5% 



neighbors. A larger n n means a greater possibility that the local velocity distribution differs 
from the overall distribution. The probability that n n is larger than the observed value K° bs , 
P{n n > K° bs ), can be estimated by Monte Carlo simulations by randomly shuffling velocities. 
Table 1 gives the results of /c-test for 215 known member galaxies, and 10 3 simulations are used 
to estimate P(n n > K° bs ) for all cases. 

It is found that the optimum scale of the nearest neighbors is 10, and the substructure 
appears most obvious on this scale. The bubble plot in Fig. 3 shows the location of localized 
variation using neighbor size n = 10, and the bubble size for each galaxy is proportional to 
— \og[PKs(D > -Dobs)] - Therefore larger bubbles indicate a greater difference between local and 
overall velocity distributions. A comparison with Fig. 2 shows that the bubble clusterings near 
the central regions of A399 and A401 are found to be significant, which indicates a distinct 
difference between the localized and whole velocity distributions. 

3 THE KMM PARTICIPATION INTO TWO CLUSTERS 

For studying the dynamical properties for each cluster, those 215 galaxies should be correctly 
partitioned. It is relatively easy to partition the galaxies whose projected locations appear close 
to the cluster centers. However, for the galaxies located exactly between the cluster centers, the 
partition might become a rather ambiguous task. 

In order to achieve a robust partition, we apply a prevalent technique of mixture modeling, 
namely KMM algorithm, on the sample of 215 galaxies. The KMM is a maximum-likelihood 
algorithm which assigns objects into groups and assesses the improvement in fitting a multi- 
group model over a single group model (Ashman et al. 1994). A detailed description of the KMM 
algorithm can be also found in Nemec & Nemec (1993). For a dynamically relaxed cluster, the 
line-of-sight velocities of galaxies are expected to be Gaussian distributed. Since A399 and 
A401 are two gravitationally distinct clusters of galaxies, we here apply the KMM algorithm 
which estimates the statistical significance of bimodality based on the three-dimension data: the 
projected positions of the galaxies and the radial velocity, just as Colless & Dunn (1996) did. 
When an initial partition into two clusters or a set of initial parameters for each cluster is given, 
the KMM algorithm can start iterating toward the maximum-likelihood solution. Table 2 lists 
the various initial partitions/parameters that we used and the corresponding final solutions, 
where {x\,y\,v\) and (x~2,y~2,V2) are the mean positions and velocities of A401 and A399, 
respectively, (a xi , <r yi , a vi ) and (a X2 , a y2 , a V2 ) are their standard deviations in positions and 
velocities, and /i and fa are the fractions of the sample in the two clusters. Furthermore, the 
estimate of the overall correction rate (R) given by the KMM algorithm is listed. 

We chose to specify initial positions and dispersions for two clusters in case 1, while in 
case 2 we specify a simple partition of sample in which all galaxies with the declination offset 
y > — 5 arcmin are assigned to be A401 members. With the different initial parameters we 
specified, the KMM algorithm searched for an optimum two-group solution, and converge to a 
very similar results. The estimate of the correct allocation rate reaches 95%. In case 3 we tried 
another initial partition: A399 concentration includes all galaxies within an angular distance of 
20 armin to the central cD galaxy UGC 2438, and we obtain the same final solution. 
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Table 2 Initial parameters and final solutions of the KMM algorithm 



Case (xi,yi,vi) 


\&X\ , O yi , O vl j 




/__ _ _ \ 
\Ox 2 ■> °V2 i ° v 2 ) 


l/i, JV 


rtate 


1 (5.9,16.9,22133) 

2 (4.0,14.6,22080) 

3 (5.0,16.8,22126) 


(9.8,9.9,1208) 
(10.7,11.1,1212) 
(10.5,9.7,1204) 


Initial Parameters 
(-11.4,-15.7,21536) 
(-11.7,-18.3,21505) 
(-9.0,-18.1,21378) 


(9.8,9.5,1227) 
(10.3,7.7,1234) 
(10.4,8.9,1232) 


(0.521,0.479) 
(0.595,0.405) 
(0.530,0.470) 




1 (4.8,14.4,22107) 

2 (4.9,14.5,22107) 

3 (4.8,14.4,22107) 


(10.2,11.7,1185) 
(10.2,11.7,1185) 
(10.2,11.7,1185) 


Final Solutions 
(-12.6,-17.4,21477) 
(-12.6,-17.3,21477) 
(-12.6,-17.4,21477) 


(9.1,8.8,1241) 
(9.2,8.9,1240) 
(9.1,8.8,1241) 


(0.586,0.414) 
(0.585,0.415) 
(0.587,0.413) 


95.2% 
95.2% 
95.1% 



According to the final partition into two clusters, there are 127 galaxies belonging to A401, 
and 88 galaxies for A399. The spatial distribution for each cluster is plotted in Fig. 4(a). Then, 
we apply the K-tests again for individual clusters, and the probability P{n n > «^ bs ) is estimated 
for each cluster. Table 3 presents the K-test results, and the corresponding bubble plot are given 
in Fig. 4(b). Compared with the last re-test on the whole binary system (see Table 1 and Fig. 3), 
the variation of localized velocity distribution for each cluster is significantly decreased, which 
indicates that the KMM algorithm arrived at a robust partition. 



Table 3 Results of re-tests for 88 galaxies in A399 and for 127 galaxies in A401 



Group size n 


2 


3 


4 


5 


6 


7 


8 


P(K n > K °„ bs ) for A399 
P(k„ > < bs ) for A401 


33.5% 
7.9% 


8.7% 
59.4% 


17.5% 
59.0% 


56.7% 
87.3% 


52.7% 
93.1% 


57.1% 
87.6% 


59.1% 
72.0% 



4 MERGER DYNAMICS BETWEEN A399 AND A401 

4.1 Velocity Structure 

The radial velocity distributions for the binary system and each cluster are given in Fig. 5. The 
solid lines represent the best-fit Gaussians with the mean velocities and dispersions listed in 
Table 2. To characterize the kinematical properties of clusters of galaxies, two robust estimators 
analogous to the velocity mean and standard deviation, namely the biweight location (Cbi) and 
scale (Sbi), are defined by Beers et al. (1990). These two quantities are resistant in the presence 
of outliers, and robust for a broad range of probable non-Gaussian underlying populations. Using 
the ROSTAT software, we compute the biweight location and scale in line-of-sight velocity 
distribution for each cluster. As a result, we obtain Cbi = 21491 ± 141 km s _1 and Sbi = 
1315±82kms- 1 for A399, and C B i = 22069±107 km s" 1 and S B i = 1212±71kms" 1 for A401. 
The physical line-of-sight velocity difference between these two clusters is V r = A(Cbi)/ (1+^) = 
539 ± 165 km s" 1 . 

4.2 Virial Mass Estimate 



A399 and A401 are gravitationally bound with each other. In order to determine the dynamical 
state of the binary system, we will apply the virial theorem for estimating the mass of each 
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cluster. Assuming that the galaxy cluster is bound and the galaxy orbits are random, the virial 
mass (M vt ) can be estimated from the following standard formula (Geller & Peebles 1973; 
Oegerle & Hill 1994): 



where cr r is the line-of-sight velocity dispersion, D is the cosmological distance of the cluster, 
A^ p = N(N — l)/2 is the number of galaxy pairs, and 6ij is the angular separation between 
galaxies i and j. The extended X-ray emission associated with A399 and A401, revealed by the 
ROSATRRl imaging observations (Fabian et al. 1997), indicates that at least the first of these 
assumptions is reasonable. We derive the virial masses of 2.0ft- -1 x 10 15 M Q and 2.1ft _1 x 1O 15 M 
for A399 and A401, respectively. 

4.3 Two-Body Models 

With the estimate of the mass for each cluster, we can investigate the state of evolution. The 
concise dynamical model for this binary system is the two-body model which was first applied 
to clusters by Beers, Geller & Huchra (1982). In this model two clusters are treated as point 
masses following a linear orbit under their mutual gravity. They are presumed to have started 
with zero separation and then to have moved apart before turning around. Given the projected 
separation of this binary system i? p , the relative radial velocity V r and the total mass M, the 
model can speculate the projection angle a (the angle between the line joining the two clusters 
and the plane of the sky), the true separation R, and relative velocity V. The equations of 
motion are as follows: 



where R m is the maximum expansion, M is the total mass of the binary system, t is the age 
of the universe, and x is the developmental angle tracing the merger process. The two clusters 
have zero separation when x — 0, 27r, while they are at maximum expansion when x — ft, 37r. 
Due to the ambiguity in observing the system only in projection, this model usually results in 
more than one orbital solution. 

The KMM analysis provides the initial estimate of the projected separation R p of the two- 
body model. Assuming a Friedmann-Robertson- Walker cosmology with fl m = 0.3, fl\ = 0.7, 
we adopt the age of the universe as to = 9.43ft -1 Gyr = 2.976 ft -1 x 10 17 s, and the angu- 
lar separation between the centroids of these two clusters (~ 36.3 arcmin) corresponds to 
a projection distance R p of 2.05ft -1 Mpc at the average redshift. We take the total mass 
M = J2 = 4.1 ft -1 x 10 15 M© in our modeling. 

Previous applications of the two-body models tried to solve the dynamical solutions within 
a range of < x < 27r, assuming that the subclusters start to saperate at t=0 and they are 
moving apart or coming together for the first time in their history (Beers, Geller & Huchra 1982; 
Oegerle & Hill 1994; Colless & Dunn 1996). However, numerical simulations of cluster collisions 
by McGlynn & Fabian (1984) showed that one clusters can even pass through each other without 
destorying the optical components. For this pair of clusters, recent observational evidences in 
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X-ray and radio bands support such a picture that A399 and A400 have been passed through 
each other in the past. As mentioned in §1, Fabian et al. (1997) found a linear X-ray structure 
in A399 pointing at the cD galaxy of A401, indicating a past violent interaction. On the other 
hand, the absence of a cooling flows in both A399 and A401 was first found by Edge et al. 
(1992), and confirmed by Markevitch et al. (1998) and Sakelliou & Ponman (2004). According 
to the numerical experiments of McGlynn & Fabian (1984), the structure of a cooling flow can 
be disrupted by the merger of two similar clusters. A simulations by Burns et al. (1997) showed 
that mergers of clusters with a mass ratio of 1:4 may destroy a pre-existing cooling flow. The 
picture of a previous interaction between these clusters is also supported by the temperature 
map for this system which shows a bridge of matter between clusters (Markevitch et al. 1998). 
Furthermore, the extensive cluster radio helo is found to be associated with A401, suggesting 
the coalescence of clusters (Harris et al. 1980). 

Based on above observations, we suppose that the two clusters started at t = with zero 
separation, and expanded to a maximum extent, and once experienced a close encounter, which 
indicates a \ range from 2it to 47r. From the equations of motion given above, four two-body 
models are allowed within 27r < \ < 4tt, including two expanding (outgoing) and two collapsing 
(incoming) models. Fig. 6 shows the (a, T4)-plane where four solutions at V T = 539 ± 165 km s -1 
are found in the bound region. Fig. 6 also shows the limit curve for the bound region, which can 
be defined by the Newtonian criterion V^R P < 2GM sin 2 a cos a. Meanwhile, we derive three 
solutions within < \ < 2n for reference only. The seven solutions are presented in Table 4. 

Table 4 Gravitationally bound solutions for the two-body model 



Case 


Dynamical 
Status 


X 

(ir rad.) 


a 

(degree) 


V 

(km s -1 ) 


R 

(ir 1 Mpc) 


(ft -1 Mpc) 


(a) 
(b) 
(c) 


Outgoing 
Incoming 
Incoming 


0.776 
1.175 
1.636 


8i.6tr 6 

75.8±i; 2 3 
9.0l 2 j 


(54CC+1B5 

556 

3460l 2 5 


14 09 +1 '° 4 

o OO+0.75 
8-.58_0.71 

2-o8±8:8_ 


15.99t_:w 
9.04l°;« 
7.101°;°? 


(d) 
(e) 
(f) 
(g) 


Outgoing 
Outgoing 
Incoming 
Incoming 


2.374 
2.854 
3.141 
3.523 


9-0 +2 2 t 
67.5±g-| 
64.4li : i 
10.3+.|i 


3430+^ 2 
583l 179 
598if 8 2 7 
3008135 


2.08l° J ;° 1 2 

5-36±8:8l 

4.74lg; 2 3 
2.081°;°? 


6.751°;°° 
5.65l°;l 2 
4.98l°;° 9 7 
4.48l°;° 



Fig. 6 shows that A399 and A401 are very likely to be gravitationally bound unless the 
projection angle a is smaller than 7°. The unbound orbit requires the true relative velocity 
V > 4, 400 km s -1 , which will lead to a very quick separation of the clusters. We should assume 
that we are not viewing the cluster system at such a special time when the projected separation 
rate reaches more than 4 ft -1 Mpc per Gyr for a pair of clusters with only 2.05 ft -1 Mpc apart. 

For cases (d) and (g), the present relative velocity of this bound system exceeds the physical 
velocity dispersion within each cluster, and the clusters are so close together that they should 
have begun to coalesce or have just coalesced. If the system were at such evolutional situations, 
we would expect to see some strong merging distortion between these two clusters in the X-ray 
surface brightness contours, contrary to the ROSAT PSPC image in Fig. 1 of Fabian et al. 
(1997). Therefore these two solutions can be definitely ruled out. 

Cases (e) and (f) are two solutions with larger projection angles, but their dynamical states 
are different. Case (e) is an expanding (outgoing) model in which the last encounter event 
occured about 2.5ft -1 Gyr ago, and the cluster pair will expand for another 1.0ft -1 Gyr, reaching 
a maximum extent of 5.65ft -1 Mpc. The true relative velocity is 583 km s -1 , and they are 
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~ 5.4/i _1 Mpc apart along the direction with a projection angle of 67°. 5. A collapsing model 
defined in case (f) is also allowed for this system. According to this model, the two clusters 
passed through each other 3.7/i _1 Gyr ago, and reached the last maximum expansion of 5.0/i _1 
Mpc about 0.8/i _1 Gyr ago. 

It is rather hard to determine whether the system is collapsing or expanding at present. The 
high resolution observations by the XMM-Newton observatory confirmed the enhancement in 
X-ray flux and temperature in the region between the two clusters, but no clues of intracluster 
compression or shock wave were found (Sakelliou & Ponman 2004). Gastaldcllo ct al. (2003) 
pointed out that the profiles of X-ray surface brightness, temperature and metallicity will shed 
light on the large-scale dynamics of the binary cluster system. Sakelliou & Ponman (2004) 
presented a clear contour plot of the residual smoothed images for these two clusters (see figure 
6 therein) where positive residuals can be found on the south-west of the central cD galaxy 
in A401 as well as on the north of the one in A399. This indicates that A401 is moving from 
south-west to the north-east while A399 is moving to the south. This feature seems to favor 
the scenario in which this pair of clusters is currently expanding. Since the projection angles 
are 67°. 5 in case (e), we can not expect to see a significant azimuthally asymmetric surface 
brightness for each cluster in the ROSAT PSPC brightness contour maps of this cluster pair 
(sec figure 2 in Markevitch et al. 1998). However, a steeper gradient can be marginally detectable 
in the north-east region of A401, and this also points toward an expanding picture. Therefore 
case (e) is the more likely solution of the two-body model. 

It should be noted that the two-body model disregards any angular momentum of the 
system, and it ignores the distribution of the matter within individual clusters which will be 
important when two clusters are close to merger. The gravitational interaction from the matter 
outside the cluster pair is also neglected. It is well appreciated that dark matter mass dominate 
in the overall dynamical mass of individual clusters. Since our estimate of dynamical mass is 
based on the virial theorem, two-body model assumes that the dark matter within a cluster is 
in approximate virial equilibrium. Despite the above-mentioned restrictions, two-body model is 
still a concise approach which is widely used to discuss the dynamic state of some gravitationally 
bound systems. 

5 SUMMARY 

We have investigated the dynamics of the cluster pair A399/A401, using the existing redshift 
measurements. We applied the KMM algorithm on a sample of 215 galaxies with known radial 
velocities, and obtained a robust separation of this cluster pair. Based on the velocity structure 
and virial mass estimate of individual clusters, we explored the two-body model for studying 
the merger dynamics between clusters. Because the observational features in the radio and X- 
ray wavebands suggest that this pair of clusters have once experienced a close encounter, we 
derived four gravitationally bound solutions within a x region between 27r and A-k. The recent 
X-ray data from the ROSAT and XMM-Newton observations can be used to choose the most 
likely two-body model. In this scenario, the pair of clusters with a true separation of 5.4/i _1 
Mpc is currently expanding at 583 km s -1 along the direction with a projection angle of 67°. 5, 
and it will reach a maximum extent of 5.65 h~ 1 Mpc in about 1.0/i _1 Gyr. 
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Fig. 1 Distribution of the radial velocities for 240 galaxies with known spectroscopic 
rcdshifts. The bin size is 1000 km s _1 
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Fig. 2 Spatial distribution for 215 member galaxies of the binary system of galaxy 
clusters, superposed by the contour map of the surface density where the smoothing 
Gaussian window with a radius of 2' is used. The contour levels are 0.03, 0.09, 0.15, 
and 0.21 arcmin -2 , respectively. 
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Fig. 3 Bubble plot showing the degree of difference between the local velocity distri- 
bution for groups of ten nearest neighbors and the overall distribution of 215 known 
cluster galaxies. 
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Fig. 4 (a) The projected positions for the member galaxies of A399 (denoted by 
astroids) and of A401 (denoted by plus sign). The dotted ellipses are the 2a contours 
of the fitted Gaussians. (b) Bubble plots for groups of six nearest neighbors for 127 
galaxies in A401 and 88 galaxies in A399. The dashed line separates two clusters. 




Fig. 5 The velocity distributions for the galaxies in (a) the entire binary system, (b) 
A401, and (c) A399, superposed with the fitted Gaussian distributions. 




Fig. 6 Projection angle a as a function of relative radial velocity V r predicted by the 
two-body model. Four bound solutions are presented at V r — 539 ± 165 km s _1 within 
2ir < x < 47r. The limit between the bound and unbound regions is also given. 



